On the Quantum Hall Effect in graphene 
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Quantum Hall efltect in 1,2-layer graphene is analyzed. The transverse and longitudinal resistivity 
are found to be universal functions of the filling factor and temperature. At fixed magnetic field 
mode the magneto-transport problem is resolved in the vicinity of the Dirac point. 
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I. INTRODUCTION 

Recently, a great deal of interest has been focused on 
the electric field effect and transport in a two-dimensional 
electron gas(2DEG) system formed in graphene flake [l|. 
In the present paper, we are are mostly concerned with 
the magnetotransport in monolayer and bilayer graphene 
in a strong magnetic field. 



II. ELECTRIC FIELD EFFECT IN GRAPHENE 
AT B=0 

We first focus on the transport properties of single- 
layer graphene (SLG) at B=0. The typical experimental 
setup is shown in Fig[Tl a. The source and drain termi- 
nals of the graphene sample are connected to the current 
source. The metal backgate and the drain are connected 
via the source of voltage, Ug, serving to change the car- 
rier density in graphene. At low-current mode / 0, the 
voltage drop, Usd-, between the source and drain contacts 
is measured, and, hence, it can be used to determine the 
graphene resistance. 

According to Refs.0,[l], the SLG energy spectrum 
obeys the linear dependence 



E±{k) ±E{k) = ±nvk, 



(1) 



where E{k) is energy, and fc, the distance in the k-space 
relative to the zone edge ( see Figll] b ). The "±" sign 
refers to the electron (+) and hole (-) conducting bands, 
respectively. The SLG state E' = is named the Dirac 
point(DP). It will be recalled that the Fermi energy, /i, in 
graphene can be varied either by means of the backgate 
voltage via field-eff'ect[H or by chemical doping. When 
/i > (fi < 0), the Fermi level falls within the electron 
( hole ) conducting band, respectively. For the special 
case, /X = 0, in which the Fermi level coincides with the 
Dirac point, the density of conducting electrons is equal 
to that of holes. 

At low currents, the voltage drop across the graphene 
sample is small as compared to applied gate voltage. Ne- 
glecting the voltage difference between the source and 
drain terminals, we plot in FigslU c-e the transverse en- 
ergy diagram for an arbitrary gate bias. At a fixed gate 
voltage, Ug, the quasi-Fermi level of the backgate contact 
is shifted by the energy, ellg, with respect to that of the 
drain contact and, respectively, to the quasi-Fermi level. 
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FIG. 1: (a)Experimental setup for gated graphene. (b)Band 
structure at k ~ 0, which shows the Dirac cones. The energy 
diagram for (c)Ug = 0, (d)[/j < and (e)f/g > 0. Xm is 
the metal work function, and Xi the electron affinity of the 
insulator. 



/i, of the graphene monolayer. Actually, the quasi-Fermi 
level of graphene sample is, pinned by that of the drain 
electrode. 

At the same time, the combined structure in FigHJ a 
can be regarded as a capacitor in which the graphene 
monolayer plays the role of one of its plate. At a fixed 
gate voltage Ug , the charge density in the graphene mono- 
layer is Q = CUg, where C = eoe/d is the capacitance 
per unit area; d, gate thickness; eo and, e. permittivity 
of free space and the relative permittivity of the SLG 
substrate, respectively. It is noteworthy that, at a cer- 
tain Fermi energy, graphene is in charge-neutrality ( CN 
) state when Q, Ug = 0. Without chemical doping, the 
state in graphene, associated with the Dirac point, coin- 
cides with the CN state. Doping shifts the DP state with 
respect to the CN state. The gate bias required to reach 
the DP state is called the offset voltage. For simplicity, 
we further neglect the chemical doping. 

Using the Gibbs statistics, we can distinguish the com- 
ponents of the thermodynamic potential for electrons, 
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FIG. 2: Conductivity of SLG at T=100;273K and fieh = 2.5 x 
lO^cm^/l/s vs. the charge density. 



f2ej and holes, flh, and, then represent them as follows 
fig = -fcT^ln h +e^^j , (2) 

k 

which gives the electron (hole) concentration as: 

Using the SLG density of states, D{E) = ^^i^i , which 
includes both the valley and spin degeneracies, we obtain 

N = 7Vo2e'Fi(l/C),P = 7Vo2e'Fi(-l/e), (4) 

where ^ = kT/ fj, is the degeneracy parameter, and Fn{z), 

the Fermi integral. Then, Nq = i {-^)^ is the density of 
degenerate electrons (holes) at T = 0. It is noteworthy 
that, at fixed /i and finite temperature, the electron (hole) 
density exceeds Nq. For example, for Dirac point state 
we obtain the intrinsic electron (hole) densities as Ni = 

f (1^) ■ '^^ instructive to estimate the typical Fermi 
energies and carrier densities. For the Fermi energy fi — 
0.2eV and v = 8x lO^m/s, we obtain the electron density 
as Nq = 4.6 X 10"'^^cm~^. The temperature range T = 
4 — 273K corresponds to a degeneracy parameter ^ = 
0.002 - 0.11. At T = lOOK, the intrinsic density is A^, = 
7 X lO^cm-2. 

We now can calculate the SLG charge density, Q, and 
the conductivity, cr, as follows: 

a = eMeh(iV + -P), Q^eiN-P). (5) 

For simplicity, we further assume a constant mobility, 
fieh, for both kinds of carriers. In Fig (21 we plot the de- 
pendence of the conductivity on the SLG charge density 



Q/e = CUg/e. As expected, the curve a{Q) is sym- 
metric and, moreover, it exhibits the minimum conduc- 
tivity, amin = '^eiJ-ehNi ^ , at the Dirac point. In 
the vicinity of the DP, when |<^| <C 1, the conductivity 

c = cTmin (1 + ^*^^): where a ^ 0.1 is a coefficient. 
For the strongly degenerate case |^| ^ 1, we obtain 

= ^ieh\Q\ ^ \Ug\. 

III. IQHE IN GRAPHENE 
A. General formalism 

We now follow the argumentation first put forward in 
Refs.jl, Q, and then modified regarding the magne- 
totransport problem of 2D electron systems in strong 
magnetic fields, including Shubnikov-de Haas Oscilla- 
tions (SdHO) and Integer Quantum Hall Effect (IQHE) 
modes. In particular, this approach allows us to repro- 
duce the SdHO beating pattern known to occur at a 
finite zero-field valley splitting (Si-MOSFET 2D system) 
and in both the crossed- and tilt-field configuration case. 

In the presence of a magnetic field, the current density, 
j, and the energy flux density, q, are given by 

j = a(E-aVr), q= (aT-C/e)j-«:Vr. (6) 

Here, E = VC/e is the electric field; C, = ^ — eip, elec- 
trochemical potential; VT, temperature gradient; and a, 
thermoelectric power. Then, a and k = LTa are the 
conductivity and thermal conductivity tensors respec- 
tively, and L = is the Lorentz number. Equation 

(|6]) is valid for a confined-topology sample, for which the 
diamagnetic surface currents Q are taken into account. 
Both the Einstein and Onsager relationships are satisfied. 

As first demonstrated in Ref.Q, the longitudinal cur- 
rent Ix = I fiowing in the circuit with an embedded Hall- 
bar sample(see Figl3]) causes a certain temperature gra- 
dient Va;T in the sample bulk. The gradient is due to 
the heating(cooling) at the first(second) sample contact, 
caused by the Peltier effect. Using Eq.(l6|), for a homoge- 
neous sample, V/i = 0, we obtain the x, y-componcnts of 
the electric field: 

Ex - aS/xT = Pxxjx, (7) 

Ey — Pyxjx • 

Here, we assume that the transverse current and tem- 
perature gradient are absent in the sample bulk, i.e. 
jyj'^yT ^ 0. Then, Pxy(x) = J^+li are the compo- 

xx~ yx 

nents of the resistivity tensor p = . 

We now calculate the longitudinal temperature gradi- 
ent V xT. The Peltier heat is known to be generated 
by a current crossing the contact between two different 
conductors. At the contact (e.g., contact "a" in FigEJ, 
the temperature T^, electrochemical potential Q, normal 
components of the total current /, and the total energy 
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flux are continuous. By contrast, there exists a differ- 
ence Aa = a^n — ot between the thermoelectric powers of 
the metal and the sample, respectively. For Aa > 0, the 
charge crossing contact " a" gains an energy e AaTf, . Con- 
sequently, Qa = /AaTa is the amount of the Peltier heat 
released per unit time in the contact "a". For Aa > 
and the current flow direction shown in FiglSJ contact " a" 
is heated and contact " b" is cooled. These contacts arc at 
different temperatures, and AT = — Tt > 0. At small 
currents, the temperature gradient is small and Ta^b ~ T . 
Therefore, the thermoelectric power of the sample, a, 
can be assumed to be constant, and we can disregard 
the Thomson heating in the sample bulk, ITS/ a = 0. It 
is noteworthy that, for conventional 2DEG, the cooling 
caused by the heat leak via the leads and the 3D substrate 
is negligible^ in quantizing magnetic fields. In fact, we 
can consider the 2D sample under adiabatic cooling con- 
ditions, when the amount of the Peltier heat released at 
contact "a" is equal to that absorbed at contact "b". Ac- 
cording to Eq.® the energy flux q is continuous at each 
contact, and, therefore, we obtain the x, y-componcnts of 
the energy flux in the form: 



(8) 



Here, we take into account that, in a finite magnetic field, 
the current is known Q to enter and leave the sample 
at two diagonally opposite high-field corners( spots ) of 
the Hall-bar sample( Figl3]). Accordingly, j'^^^^ are the 
X — y components of the current density in the vicin- 
ity of the high-field corner. Evidently, the total cur- 
rent that enters(leaves) the sample at high-field spots. 



' {I'xY + i-^yYi equal to the current in the sam- 
ple bulk, i.e. I' = I = j.j;d, where d is the sample width. 
Folding Eq.®, we finally obtain the longitudinal tem- 
perature gradient VrT = ^^^=, which is linear 

in current. 

Using Eq.©, the voltage drop U, measured between 
the open ends "e" and "d" can be easily found a.s U — 
J Exdx = Aa{Ta — Tf,) + Pxxjxl, where I is the sample 
length. With the resistance and the thermoelectric power 
of the metal leads disregarded, the total resistivity p = 
U/jxl and the Hall resistivity pn ~ Ey/jx of the sample 
are given by 




FIG. 3: Experimental setup. The red arrows show the current 
flowing at the high-field corner and in the sample bulk. 



exceeds both the thermal energy and the energy re- 
lated to the LL-width. 2D carriers can be considered 
dissipationless [lo| in this case, and, therefore, we assume 
<^xxtPxx — and thus 9h ~ tt/2. Equation ([9]) yields the 
result 



PH 



(10) 



previously reportedQ for the unipolar magneto- 
transport by 2D electrons or holes. 

In contrast to the conventional 2D systems, carriers 
of both kinds, electrons and holes, coexist in graphene 
[H [ll|. Therefore, Eq. lfTO)) is only valid when the 
current-carrying state is provided by the majority car- 
riers (electrons N ^ P or holes N P). We emphasize 
that in strong magnetic fields the Hall conductivity of 
the electron-hole plasma 



{N - P)ec 
B ' 



(11) 



may be low enough and, moreover, it vanishes at the 
Dirac point when N = P. The Hall angle vanishes as 
well, 9h ^ 0- Actually, the current flow at DP is uniform 
throughout the sample, including the contacts. The high- 
flcld corners disappear. In the ultra-narrow vicinity of 
the DP, i.e., at ^ a^x Eq.® yields 



p^aJil + a^L), 



PH 



(12) 



P = Pxx + 



a 



(JyxL 



sm Wff, 



PH = Pyx, 



(9) 



where Oh = arctan(crj,a;/(Ta;^) is the conventional Hall an- 
gle. It is noteworthy that the total resistivity is an even 
function of the magnetic field. Equation © is the cen- 
tral result of the present paper. First, at zero magnetic 
field, we obtain Oyx, Oh = 0, and, therefore, Eq.® repro- 
duces the result p = cr^^(l + /L) reported by Kirby 
and Laubitz Q. 

The primary concern of the present paper is the case 
of strong magnetic field, in which the cyclotron energy 



According to Ea. (|T2)) the total resistivity p is finite, 
whereas the Hall resistivity vanishes in agreement with 
experiments [l| . 



B. IQHE in SLG at variable magnetic field 

In a strong magnetic field, the SLG spectrum is given 
byi: 



hv 



(13) 
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where Xb = \/hc/eB is the magnetic length, and = 
0, 1,2.., the Landau level(LL) number. Here, the ± sign 
refers to the electron(+) and hole(-) energy bands, re- 
spectively. We assume further no LL broadening. 

With the energy spectrum specified by Eg. (fT3| . the 
electron(e) and hole(h) part of the thermodynamic po- 
tential fle.h are given by 



-fcrr 



ln(l-t- ) -I- 2 ln(l -I- ) 



JV=1 



,(14) 



where F = (27rA^) ^ is the zero- width LL density of 



states. In Eq. (|T4| . we take into account the four- fold 
degeneracy( spin-|-valley ) for high-index {Af > 1) LLs. 
Then, the zeroth LL exhibits Q a double degeneracy for 
both electrons and holes. 

Let us introduce the dimensionless electron energy 
spectrum e_\f = Ej^/fj, = sign{i') 4J\f where v = 
sign(^)A'o/r is the conventional filling factor. We now 
consider as a starting point the magneto-transport solely 
by electrons {N ^ P) at a finite temperature and cer- 
tain Fermi energy /i = Hvy/nNo > 0. Using the thermo- 
dynamic definition specified by Eq.® and Eq. p^ . we 
finally obtain the electron density 



= F 



2T 



(15) 



where T is the Fermi function. 

It will be recalled that, in a str ong magnetic field, the 
thermoelectric power is known Il2| to be proportional 
to the entropy 5*6 per particle. For electrons, we obtain 
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FIG. 4: Magnetic field dependence of the total and Hall re- 
sistivities, calculated using the energy spectrum Eq. (|13p and 
Eas. (|10ll5ll6p for electrons(/i > 0) at ^ = 0.005,0.01,0.02. 
The thin line represents the low-B asymptote pH ~ 
The dashed line represents the universal resistivity asymp- 
tote p""™. Inset: enlarged plot of the low-B p-data 



power, carrier density and, hence, resistivities given by 
Eg. pT]]) approach the universal values ai"^' = 



pf^ = ^^J^\ andp(-^) = pf\af'>Y/L, irrespective of 
temperature. The dotted line in Fig |4] represents the uni- 
versal resistivity asymptote p""™(i/) = ^ ^''^J^F'Ji whose 
values aX V = vj^ correspond to the universal resistivity 
values Then, using Eqs. (|10|16l27p . we plot in the 

inset of FigUlthe low-B magnetorcsistivity which exhibits 
onset of SdH oscillations at ~ ^. We emphasize that 
the SLG energy spectrum specified by Eq. (|T3| suggests 
that zcroth-LL always remains occupied at /i > 0, and, 
therefore, the Hall resistivity is finite even in ultrahigh 
magnetic fields. 



It is noteworthy that the thermodynamic potential ilg, 
density A^, and thermoelectric power are universal 
functions of the dimensionless magnetic field and 
temperature ^. In the case of low temperatures and weak 
magnetic fields z^^^,^ ^ 1, we use Lifshitz-Kosevich for- 
malism and, thereby, derive in Appendix |VA] asymptotic 
formulas for 57e, N ^ S'e, and, hence a^^. 

Using Eq. (|10|15ll6p . in Fig|4]we plot the dependencies 
of the Hall resistivity = and the total resistiv- 
ity p as a functions of the magnetic field. At low mag- 
netic fields, ph ~ ^ B. Then, in a quantizing mag- 
netic field, the Hall resistivity exhibits the QH plateau 
at PH = 4j)+2 ' accompanied by a plateau-to-plateau 
transition at the so-called critical filings, vj^f = AN, 
when the Fermi energy coincides with the A/'-th Lan- 
dau level. Simultaneously, in the vicinity of the criti- 
cal fillings vj\f, the total resistivity p exhibits peaks. It is 
noteworthy that, at critical filings i^a/, the thermoelectric 



C. IQHE in SLG at fixed magnetic field 

We now examine the magneto-transport in SLG at a 
fixed magnetic field for an arbitrary /i, i.e., in presence of 
electrons and holes. At a moment, we neglect the longitu- 
dinal conductivity of the 2D electrons i.e. Oxx = 0. Using 
Eqs. (|3|14|) and, in addition, combining the electron(hole) 
carrier contributions, we finally obtain the Hall conduc- 
tivity of monolayer graphene 
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(17) 



According to Eq. (111117^ , the Hall conductivity is an odd 
function of the magnetic field and Fermi energy. In FiglSj 
the Hall conductivity is represented as a function of the 
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FIG. 5: Hall conductivity dyx and the total resistivity p 
given by Eq.(Il7| and Eqs. (|10I17I18P . respectively, plotted 
as functions of the filling factor v at different temperatures 
i?i = 0.02,0.1 for energy spectrum specified by Ea. (|13p . The 
dashed line corresponds to that in FigU In the vicinity of the 
Dirac point (shaded area), the total resistivity is plotted in 
Fig[6]for the split-off zeroth-LL energy spectrum. Inset: SLG 
thermoelectric power given by Ea. H18|) vs. the filling factor v 
at i9i = 0.06,0.1. The dashed line represents the asymptote 



where a/i is the thermoelectric power of holes. The ther- 
moelectric power of monolayer graphene, specified by 
Eq. ([T8|) . coincides with that for electrons and holes at 
N ^ P and N <^ P, respectively. For strong magnetic 
fields, we use Eas. (|10ll7ll8p and, plot in Figl5]the total 
resistivity against the filling factor. In agreement with 
experiment [H IT3| . the dependence ^(i^) exhibits peaks in 
the vicinity of the critical filings z/jv, except for zcroth- 
LL-related peak at ;^ = 0. Below, we attempt to explain 
this discrepancy. 

We argue that the magnetoresistivity vanishes near the 
Dirac point (see FiglS]) which is caused by the 4-fold de- 
generacy of the zeroth-LL filled by both electrons and 
holes. Indeed, only the zeroth-LL can make the main 
contribution to the electron il^ = 2rin(l -I- ) and hole 
rih = f2e(— thermodynamic potentials. After some al- 
gebra based on Eas. (|18l3p . we can surprisingly find that 
a = (see the inset of FiglS]) within a wide range of Fermi 
energies when | |< 4. Accordingly, the total resistivity 
vanishes within the above i^-rangc as well (see the shaded 
area in Figl5|). 

Our results seem to contradict to the experimental 
observations [H ITsj . according to which a well-pronounced 
p-peak is observed at the Dirac point. We now demon- 
strate that the above discrepancy may be due to lifting 
of the four-fold degeneracy of zcroth-LL in real systems. 



filling factor v. It is noteworthy that the dimensionless 
temperature embedded into Ea. ((T7| depends on the fill- 
ing factor. Indeed, ^ — — = _ , where we in- 

M sign(iy)^|i/| 

troduced the dimensionless temperature at a fixed mag- 
netic field, i?i = ^^^^^ ■ For the typical B = 9T, 
V = 8x lO^cm/s, and T = 4.2K, we estimate = 0.01. 
Figure [S] demonstrates the sequence of the QH plateaus 
ayx — ±^{J\f + 1/2) accompanied by the plateau-to- 
plateau transitions at the critical values a^'^ = ^VN' 
of the Hall conductivity. In the vicinity of a certain 
plateau-to-plateau transition, /S.v ~ v — vj^^^v <C vj^^ 
the Hall conductivity Oy^ = cr^'^^''(l + ^ f)!^^ ) exhibits 

the linear T-scaling. It is noteworthy that, in contrast to 
conventional 2D systems, the width of the A/'-th plateau- 



3/2 



AT, 



3/2 



to-plateau transition in SLG increases as Vj^ '~ j»o 

Of particular interest is the behavior of the Hall con- 
ductivity in the vicinity of Dirac point. According to 
energy spectrum given by Eq. p^ . the electron ( hole ) 
zeroth-LLs are half-filled at ^ = 0. Thus, the Hall con- 
ductivity specified by Eg. lfTT)) exhibits a 4e^//i step at 
DP within the ultra-narrow range of fillings v ^ 'd\. 

In order to complete our study of the magneto- 
transport in SLG for mixed state containing both elec- 
trons and holes, we find the thermoelectric power 



a 



a^N + ahP 
N + P ' 
-ae(-/i) 



(18) 



ah 



D. IQHE in SLG in the vicinity of the Dirac point 

In this section we resolve the magnetotransport prob- 
lem in the vicinity of the Dirac point taking into account 
the zero-LL degeneracy lifting, and moreover assuming 
finite SLG conductivity a^x 7^ 0. The nature of the en- 
ergy states in graphene at A/" = LL has recently been 
under intense theoretical and experimental [l3l - l2l| inves- 
tigation, but still remains unresolved. Conventionally, 
the four-fold degeneracy of the zcroth-LL is associated 
with a two-fold degeneracy resulting from the spin sym- 
metry and a two-fold degeneracy caused by sub-lattice 
symmetry. We now make an attempt to distinguish the 
hierarchy of the split-off zeroth-LL on the basis of the 
available experimental data( see Table |I| . 

For moderately disordered samples at a fixed temper- 
ature and a fixed magnetic field < SOT, the density de- 
pendence of the magnetoresistivity exhibits a symmetric 
peak centered at DP. The height of the resistivity peak 
increases as temperature is lowered. Then, the set of 
the density-dependent resistivity isotherm curves (typi- 
cally 0.5 — 50K) clearly demonstrates two T-independcnt 
points p ^ 0.5/i/e^ on the electron and hole sides of the 
resistivity peak. These points lie at certain critical fillings 
^ ±0.5( see open symbols in FiglHlinset "a"), which we 
attribute to crossing of the electron(+) and hole(— ) sub- 
branches of the split-off zeroth-LL by the Fermi level. It 
is noteworthy that, usually, experimental /5-data are con- 
verted to give the inverse resistivity tensor p~^. At the 
above-mentioned critical fillings of ^ ±0.5, the longitu- 
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TABLE I: The device parameters used to evaluate the zeroth- 
LL splitting and the magnetoresistivity at the critical fillings 
t'o^i and ^^-2 in Fig|6]inset "a" 



B,T T,K 






4 Ref. 




8 2.3 


2 


0.6 


M 


14 0.3-100 


1.3 


0.98 


[16] 


18 0.5-50 


0.6 


0.5 


[18] 


20 4.2 


2 


0.27/0.31 


[20] 


30 4 




0.42 


[15] 


0.4-150 


1 


0.38 


[19] 


45 4.2-294 


2 


-/0.3 


[20] 


1.4 


5 


-/0.2 


[21J 


9-45 0.03 


5 




[2^ 



dinal component of the inverse resistivity tensor (p ^)xx 
also exhibits T-independent points, whereas the trans- 
verse component {p~^)xy experiences a smooth (0 — 2)j- 
transition [l3 - [l6l fisl. [l9j . Actually, the /5~^-data provide 
an alternative method for finding the critical fillings. 

Further progress related to the zeroth-LL splitting 
problem has been made for high-mobility samples at 
low temperatures [2l| . Surprisingly, at ultra- low temper- 
atures 30mK) the previously reportedflB, [l^ smooth 
transition (0 — 2)^ in {p~^)xy turns into a sequence( 
(0 — 1)^ and (1 — 2)^) of two sharp transitions. These 
transitions occur at certain critical fillings Vq_i and 
vf_2 associated with the crossing of the electron(-t-) and 
hole(— ) related lowest(O-l) and highest(l-2) zeroth-LL 
split subbranches by Fermi level. In FiglHlinset "a" we 
plot a fun diagram of split-off zeroth-LL critical fillings 
and, moreover, add the data of [13, [HI for ultra-high 
fields < 45T at which the resistivity exhibits well de- 
fined split-off zeroth-LL peaks near the DP. For the lat- 
ter high-field case, the critical fillings were determined by 
careful analysis of the split-off resistivity peak positions. 
However, exact definition of the zeroth-LL subbranches( 
spin and(or) valley ) remains unclear. The available tilt- 
field experimental dataflO, [HI point to only the spin- 
independent origin of zeroth-LL splitting. 

Of particular interest is the behavior of the density- 
dependent resistivity p measured in the vicinity of the 
Dirac point at a fixed temperature for different ap- 
plied magnetic fields. When the magnetic field grows 
[B < 31T), both the height and width of the DP resis- 
tivity peak increase [l7|. Moreover, at the Dirac point 
p{B) dependence exhibits [TgI. [TtI the abrupt rise at a 
certain critical magnetic field. The higher the sample 
mobility at i? = 0( i.e. the smaller the offset gate bias 
), the lower the critical magnetic field. Intriguingly, the 
expected divergence of the DP magnetoresistivity is not 
observed [2^ for high-disorder samples ~ lO'^Spj up to 
B = 57T. We argue that these effects are highly expected 
within our scenario. It will be recalled that the resistiv- 
ity is a universal function of the filling factor v and the 
dimensionless temperature i9i ^ kT/^/B. Consequently, 
the observed [l6| enhancement of the resistivity, caused 




FIG. 6: Evolution of the total resistivity p in the presence of 
zeroth-LL splitting at critical fillings = ±0.5 and u^_2 = 
±0.75 ( corresponds to B=14T in inset a ) for dimensionless 
temperatures i9i = 0.3,0.25,0.2,0.15,0.1,0.05,0.02,0.015. 
SLG conductivity a^x ~ O.Ole^/h. Inset a: Hierarchy of 
the split-off subbranches of zeroth-LL. Direct observation of 
critical filling, based on T-independent points: >-[3, <l-[l^. 
O-fi^- Indirect analysis of the critical filling on the basis of 
p^^ data: ■A'- [3, '-[IH. Proximate analysis based 

on the position of the split-off p-peaks A-[2^, '-[111. Inset b: 
HaU conductivity at i?i = 0.1,0.015. 



by the increase in the magnetic field at r=const, may be 
similar to that reported |17H13 | for the T — > sweep at 
i3=const. 

To compare our results with the experiment p^. [l7| . we 
plot in Figinithe transport coefficients p, cjyx in the vicin- 
ity of the Dirac point on the assumption that the crit- 
ical fillings have values i^q_i — ±0.5 and v^_2 = ±0.75 
deduced from FiglHlinset "a" for B=14T. Then, we set 
the constant longitudinal conductivity axx ~ O.Ole^/h, 
which provides a reasonable value for the sample resis- 
tivity axx ~ 2.5Mn at the Dirac point. It is note- 
worthy that, at r < 5K the p(i^)-curves found to be 
indistinguishable in ^isf , thus we can roughly estimate 
the typical LL width as ~ 5K. Finally, we set in FiglS] 
■di = 0.015-0.3 for the temperature range T = 10-225K. 
At high temperatures, the T-independent point at the 
electron(hole) sides is clearly seen. By contrast, at 
T ^ 0, the electron(hole) zeroth-LL subbranches be- 
come resolved. The ultranarrow peak(deep) at the 
Dirac point is attributed to 2D graphene resistivity when 
(Tyx ^ cTxx- We emphasize that dependencies of the 
p{B),p{T) at the Dirac point 1/ ~ can provide the 
unique knowledge regarding the own magnetoresistivity 
Pxx of 2D graphene fiake. 

It is instructive to compare the resistivities at the 
critical fillings z^(o-i) , 1^(1-2) with those reported in ex- 
periment( see Table |T] ). For half-filled lower(O-l) and 
highcr(l-2) zeroth-LL subbranch wc obtain the electron 
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densities 1/2V and 3/2r, respectively. Then, the en- 
tropy = fcrin2 is the same in both cases. The 

critical resistivities arc n^'^^^^ ~ 1' = 1.17-% and 

p*^^"^^ = 0.043^, respectively. According to Fig|6l only 
the critical point at persists at elevated temper- 

atures, thus, the respective critical resistivity p^"^^^ is 
consistent with that 0.5-%) observed in experiment 
(see Table HI). 

F inally , we calculate the absolute value of ener gy E = 
hvy/irVv or in dimension units i?(meV) = 14, h\/vB(T) 
related to zero LL split sublevels. In Figl7]we re-plot the 
results represented in Fig|6l inset "a" in terms of zero LL 
split energy vs magnetic field. 



E. IQHE in bilayer graphene 

We now analyze the case of bilayer graphene(BLG) 
magnetotransport, which differs insignificantly from that 
considered above for monolayer graphene. With our pre- 
vious notations, the BLG energy spectrum [23| is given 



by 



= ±Em = ±nw,^N{N -I), 



(19) 



where Wc = is the cyclotron frequency; m = = 
0.054me, effective mass; 7 = 0.39eV, interlayer coupling 
constant. For the general case of an arbitrary interlayer 
coupling, the BLG spectrum is described[26j by Eg. ([50)) 
(see Appendix IV B[) . 

The electron- and hole-related parts of the BLG ther- 
modynamic potential are 



-kTT 



4 ^ ln(l + e^TT^) + 41n(l + ef^) 



N=2 



,(20) 



where we take into account the four-fold degeneracy( 
spin+valley ) for M > 2-index LLs. Then, based on the 
energy spectrum given by Eq. (jl9|) we take into account 
the total 8-fold degeneracy of the zero-energy LL due to 
the spin, valley, and A/" = 0, 1 orbital LL degeneracies. 

Similarly to the SLG case, we introduce the dimension- 
less electron energy spectrum as ej\f = ^\/N[N — 1). 
Using Eqs.(Iini I3ll6ll0p we plot in FigE] the magnetic 
field dependence of the Hall resistivity pn = and 
of the total resistivity p. The dashed line represents 
the resistivity asymptote p'^'^™{v) whose values at vj^ = 
4VAA(A/'- 1),AA > 2 gives the respective universal val- 
ues when the Fermi energy coincides with the A/'-th Lan- 
dau level. 

We now extend our previous results for the bipolar 
magneto-transport in SLG at i?=const for the present 
case of bilayer graphene. Using Eas. (|20l3llip . we obtain 
for the BLG Hall conductivity 



4e^ 
~h 



E 



th 



+ 



(21) 



It is worthwhile to mention that the dimensionless tem- 
perature embedded into Eq.([2T|) can be re- written as 
£ = — , where we introduced the dimensionless tem- 
perature at a fixed magnetic field, i?2 = j^- In FiglHl 
the Hall conductivity is represented as a function of the 
filling factor i>. Assuming that the charge transport 
is dissipationless in strong magnetic fields, we employ 
Eas. (|10l3ll8l21|) and plot the total resistivity in Figlll It 
exhibits peaks in the vicinity of the critical filings i/_\f in 
agreement with experimental findings p^]. 

In order to analyze the p-behavior in the vicinity 
of the Dirac point, we use the available experimen- 
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FIG. 9; QHE in BLG at fixed magnetic field B = 20T, di- 
mensionless temperature 1)2 ~ 0.5,0.2 and a^x ~ O.Ole^/h. 
The dashed line correspond to universal resistivity asymptote 
p""™. Inset: magnetic filed dependence of the electron(+) 
split-off — 0,1 LL subbranch fillings. The arrows in the 
main panel corresponds to fillings at B = 20T in inset. 
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FIG. 10: BLG hierarchy of the split-off electron subbranches 
of TV = 0, 1 LL. The critical fillings are extracted from p~^- 
data □-[23|; analyzing the position of the split-off p-peaks 



V. APPENDIX 



Lifshitz-Kosevich formalism 



tal data[25[ and identify the split-off subbranches of 
the zero-energy LLs{J\f — 0,1). Our result is repre- 
sented in Figiniinset for electron branch(niarkcd as "-f") 
of the BLG spectrum. In strong magnetic fields the 
electron-related zero-energy LLs splits into four sub- 
levels ^Q-ii^i-2i^2-ii^t-4: shown in FiglQlinset. Simi- 
lar to monolayer graphene case the resistivity peak at 
the Dirac point is caused by split-off subbranches of the 
zero-energy LLs. 

Finally, we calculate the absolute value of energy E = 
^^^^^Tv or in dimensional units £;(meV) = 0.53z^B(T) 
related to zero-energy split LL sublevels. For electrons 
the result is shown in FigHU] As expected, the A/" = 0, 1 
LL splitting disappears at zero magnetic field. 



Using the conventional Poisson formulae 

00 00 
00 00 

^^(A/-)= / (^(AA)dAA+2Re^ / ^{N)e^^'''^ dN , 

.(22) 

where mg — 1 < a < toq ("^o is lower limit of summation) , 
the thermodynamic potential can be represented as the 
sum Qe = + of the zero-field and oscillating parts 
as 



'A^oK'J^2 (1/0(23) 



= -8 



^2jfixReX: /e2-^^^ln(l + e^)dAA, 



After a simple integration by parts, the oscillating term 
has the form 



IV. CONCLUSIONS 



In conclusion, we studied the Quantum Hall effect in 
mono- and bi-layer graphene. The transverse and lon- 
gitudinal resistivities were found as a universal function 
of the filling factor and temperature. Based on exper- 
imental data, we found the B-dependence of zeroth-LL 
split subbranches and, in addition, resolved the magneto- 
transport problem in the vicinity of the Dirac point. 



0„ 



^ fe=l ^''^'^ i 1+65 



(24) 



It should be noted that the integrand in Eq. ((24| is 
a rapidly oscillating function, which is, in addition, 
strongly damped at e > 1. A major part of the integral 
results from the energy range close to the Fermi energy, 
when e ^ \. Therefore, M{e) = ^ can be regarded as 
smooth functions of energy, and, hence, can be re- written 
as 



AAUi=- + ^(.-l). 



4 2 



(25) 
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Under the above assumption, we can change the lower 
hniit of integration to — oo and then use the textbook 

expression / jj-^dy = sinhUk\ fo'^ ^^le integral of the 



above type, 
tential 



sinh(7r/c) 
-oo 

Finally, we have for the thermodynamic po- 



E 

k=l 



k sinh ri 



(26) 



where = Tr^fcz/^ is a dimensionlcss parameter related 
to T-damping of the oscillations ^ cos(7rfcj//2). Using 
the conventional thermodynamic definition and Ea.(|24p. 
we find both the entropy and the density of Dirac elec- 
trons, valid at low temperatures and weak magnetic fields 



N = Nt+ ANqttS, ^ 



sm 



^) -h7rCcos(^)coth(rfc) 



k=l 



sinh(rA;) 
Sp = S - AknNn > ' cos 



(27) 



oo 



kv 



where Nt 



7Vo2C^i^i(l/C) and St 

are the carrier 

density and entropy at _B = 0, respectively; and, 
*&(z) = is the form-factor. 



B. 



sinh(2:) 



Spectrum of bilayer graphene at finite 
interlayer coupling 



At B=0, the spectrum of bilayer graphene spectrum 
was found in Ref.[2? ~- 



as 



E{k) 



±-7 

2 ' 



(28) 



For clarity, we further consider the electron-related 
part("-l~") of the spectrum. For a strong inter-layer cou- 
pling < 1, the BLG spectrum E{k) = ^ coin- 
cides with that known for regular 2DEG systems. In the 
opposite case of a weak interlayer coupling, i.e., a high 



carrier momentum 



2vhk 



» 1, Eq. ((28)) reproduces the lin- 



ear Dirac spectrum E(k) = hvk. For an arbitrary 7, we 
can find the BLG density of states D(E) 



2-E+7 

then calculate the electron density for the actual T 
case as 



and 



TV = iVo 1 + 



(29) 



For strong(weak) inter-layer coupling, Eq. (|29| yields the 
relationship for the Fermi energy, N = = 
known for the conventional 2DEG and that ( N = Nq) 
for Dirac electrons, respectively. The transition from the 
Fermi-like spectrum to that of the Dirac type occurs at 
22i^ ^ 1 , which yields the critical electron density N* ~ 




TIN 



FIG. 11: Phase diagram of the BLG electron states for the 
lowest(A/' = 1,2,3) LLs and an arbitrary inter-layer cou- 
pling. The red and blue lines are related to the BLG(!^ = 
4:^J\f[J\f — 1)) and SLG(j^ = AM) energy states, respectively. 
The green line represents the BLG-to-SLG energy phase sep- 
arator. 



^^fi2^j2 ■ It is noteworthy that Eg. ([29]) makes it possible to 
find explicitly the Fermi energy ^{N) = j{^/l + N/N* — 
l)/2. The BLG-to-SLG transition of the energy spectrum 
occurs at iV = N* , which corresponds to the critical 
Fermi energy fi* — 7(\/2 — l)/2. Assuming that 7 ~ 
0.39cV and u = 8 X 10^ cm/s, we can obtain the following 
estimates N* = 4.36 x lO^^ cm-2[23| and ^t* = 0.08eV. 

We emphasize that, in a strong magnetic field, the 
inter-layer coupling may play an important role as well. 
For arbitrary 7, the LL spectrum of the bilayer graphene 
was derived in Ref. 12611 
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f3-l + 2M- V(/3- l)2 + 4AA/3 



1/2 



(30) 

Here, we introduce the inter-layer coupling parameter 
/3 = iV*/r. For the uhimate case of SLG(/3 = 0) and 
BLG(/3 — !• 00), it is possible to easily reconstruct the dis- 
crete LL energy spectra specified by Eq. (fT3l) and Eg. ([T9)) . 
respectively. The BLG-to-SLG change in the energy 
spectrum occurs at /3 ^ 2{2Af —1). The higher the LL 
index, the lower the magnetic field at which the BLG-to- 
SLG transition occurs. 

We can now find the phase diagram of BLG elec- 
tron states N{T) for an arbitrary inter-layer coupling. 
Using the condition £V = /i(Af) related to the Fermi 
level crossing the Af-th LL and Eas. (|29l30p . we calcu- 
late the phase diagram( see FigfTTj). Since both axes 
are scaled in density units N*. the phase diagram is 
valid for an arbitrary inter-layer coupling 7. With in- 
creasing magnetic field strength, the initial BLG energy 
state( with a certain LL index ) is transformed to that 
belonging to SLG. We ca n easily fin d the BLG-SLG 
phase separator N/N* = y^l — (2r/A^*)^ represented by 
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the green line in Fig[TT] As expected, at B 0, the 
phase diagram demonstrates a BLG-to-SLG transition at 
N = N* . It will be recalled that, in a real BLG system, 
N* ~ 4.36 X 10^^ cm~^, and, therefore, the magnetic field 



scale of the BLG-to-SLG transition (i.e. F ^ N* /2) is on 
the order of i? ^ 90T. We conclude that, for the actual 
experimental range B ~ Q~ 40T the BLG spectrum can 
be well described by the simplified Eq. p^ . 
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